arXiv: 1509.03445vl [math.AP] 11 Sep 2015 


MOTION OF THE GINZBURG-LANDAU VORTIGES FOR THE 
MIXED FLOW WITH CONVEGTIVE FORCING 

OLGA CHUGREEVA 


Abstract. We consider the mixed Ginzburg-Landau flow that is supplemented with 
convective derivatives of the unknown function. We show that the associated vortex 
motion law is the mixed flow of the renormalized energy with new nonlinear forcing 
terms. These terms are uniquely determined by the extra terms in the initial PDE. 
Our proof relies on the assumption that the initial data are close to optimal. 


1. Introduction 

We study the mixed Ginzburg-Landau flow that contains additional convective terms. 
The equation reads 

( 1 ) (Ae -I- i)dtUe + ke{F ■ V)Ue -L (G • V){iUe) = AUg -I- 2(1 _ 

The function Ug is complex-valued. It is defined on a smooth, bounded, simply connected 
domain D C The vector fields F{x, t) and G{x, t) are smooth in both variables. The 
convective derivative {F ■ 'S/)ui; is a two-dimensional vector field with the components 
{F,Vui) for j = 1,2. The parameter £ is positive. The scaling factors are = 1 /| loge| 
and Ae = Aq • with Aq positive. We complement © with either Dirichlet or zero 
Neumann boundary conditions. 

We are interested in the asymptotic behavior of the solutions to o when e goes to 
zero. The key phenomenon is the emergence of point singularities of the solution, the 
so-called vortices. The behavior of is effectively reduced to the motion of the vortices. 
Our main result is the identification of the vortex equation that arises from ([T}. We 
show that the vortices obey the system of ODEs 

(2) (Ao -L i)ak - F{ak{t),t) - dkiG{ak{t),t) = --da^W{a{t)). 

TT 

Here, ak is the position of the fcth vortex and dk is the corresponding degree - the 
topological charge carried by the vortex. The renormalized energy W (a) depends on 
the overall vortex configuration a and encodes the interaction between the vortices. We 
note that the external vector helds F and G are transferred from the initial PDE into 
the effective motion law. 

The equation o is the natural generalization of the evolution equations based on 
the Ginzburg-Landau energy functional 

(3) E^{u) := j ^\Vu\^ + dx. 

D 

The expression on the right-hand side of © is the negative gradient of Es{ue). The 
functional Es{u) is the reduced version of the energy functional that is used to describe 
superconductivity and superfluidity. In this regard, the vortices of provide a toy 
model of the Abrikosov vortices that are observed in type-II superconductors. 

The vortex dynamics associated to various flows of the Ginzburg-Landau energy has 
been studied since early 1990es. It is well understood by now. The vortex motion law 
for the Ginzburg-Landau-Schrbdinger equation (Hamiltonian flow of Egiu^)) has been 
first derived by Neu m- The first result of this kind for the heat (gradient) flow of 
Es(us) is due to E [5]. Both authors used the method of matched asymptotics. More 
rigorous proofs for the heat flow have been given by Lin US], Jerrard-Soner [7] and 
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Sandier-Serfaty m- The Schrodinger flow has been further investigated by Colliander- 
Jerrard [5] and Lin-Xin m- The equation on vortices for the mixed Ginzburg-Landau 
flow, i.e., for o with F — G = 0, has been established by Kurzke-Melcher-Moser-Spirn 
m and by Miot [19]. In all the cases, the vortex dynamics is governed by the flow of 
the renormalized energy W. Depending on the original PDE, the flow of W is again 
gradient, Hamiltonian or mixed. 

When adding non-variational terms to a Ginzburg-Landau equation, we are interested 
in perturbations that modify the vortex dynamics without destroying it. Suppose that 
we start with a forced flow of E^iu^). It should produce an equation on vortices that 
itself is a forced flow of W. We show that the forcing with convective derivatives satisfies 
these requirements. The requirements are further motivated by the prospect of studying 
stochastic Ginzburg-Landau equations. In the stochastic setting, we are concerned with 
the effect of a random external field on the vortex dynamics. First results in this 
direction are presented in our work | 2 ]. 

The coefficient kg in o singles out the correct time scale for the vortex motion 
and the correct amplitude for the external helds. The slow motion of the vortices is 
characteristic for the Ginzburg-Landau heat flow. The Ginzburg-Landau-Schrodinger 
equation does not need the rescaling. We have to damp the F term, but not the G 
term, to see the effect of both fields on the vortex motion. This hints that the term 
kg{F ■ V)Me is related to the damped part of the time-derivative, X^dtUe, and thus to 
the heat flow. In the same way, the term (G • V){iUe) is related to the Schrodinger flow. 
If we replace in dm the usual time-derivatives with the material derivatives, we arrive 
at the equation 

XeDFU-\-iDGU = Aue + ^(1 - \Ue\^)Ue. 

For an arbitrary vector field H, Dhu := dtu + {H • V)u is the material derivative 
defined by H. From this point of view, m is again a combination of the heat and the 
Schrodinger flow. The scaling on the left-hand side of © simply reproduces the scaling 
in the mixed Ginzburg-Landau flow without forcing. This leads to the richest possible 
vortex dynamics. The equation @ takes into account both the interaction between the 
vortices and the impact of the external fields F and G. If the F term in m is damped 
stronger than by kg, it does not affect the vortex motion. Similarly, if the G term is 
damped in any way at all, it does not appear in ([ 2 ]). 

Another reason to study m is the formal analogy between the mixed Ginzburg- 
Landau flow and the Landau-Lifshitz-Gilbert equation from the theory of ferromag¬ 
netism. The closest counterpart of o is a particular variant of the latter, the Landau- 
Lifshitz-Slonczewski equation. It describes the evolution of the magnetization direction 
rUg : D —>■ in the presence of an external spin-current v. In this context, D represents 
a thin magnetic sample. The equation reads 

(4) dtirig -I- {v ■ y)mg + ^rrig x (u • V)mg = rrig x {XgdtrUg + VL 2 Eg{mg)). 

The basic mathematical properties of (Hj) have been scrutinized in [18] . 

In certain situations, the magnetic energy Eg{mg) is modeled by a functional similar 
to the Ginzburg-Landau energy. The functional is given by 

Eg{m) = J i|Vmp -f ^rnl dx. 

D 

Here, m 3 is the third component of the vector m. For a more detailed discussion, see 
Section 7 in (T^j. 

We can formally link equation (|3|) to o in the following way. We project (Hj) on the 
tangent plane that is identified with C. Then, the vector multiplication by rrtg 

becomes the multiplication by i. This transforms equation into equation o with 
F — G = V and kg = [3. 

The solutions to (|3|) can develop point singularities similar to the Ginzburg-Landau 
vortices. These magnetic vortices move according to the Thiele equation, which is 
analogous to For the case = 0, the equation has been known since the works of 
Thiele and Huber [Bj. It has been rigorously derived in [T2] and [13]. The passage 
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to the vortex motion law for (|3]) has been justified in [TU]. The role of the field v is 
essentially the same as that of the fields F and G in ([2|) . Our conclusions are therefore 
the Ginzburg-Landau counterpart of the results in m- 

The study of the Landau-Lifshitz-Slonczewski equation is, however, more delicate 
than that of the Ginzburg-Landau equation. The main problem is that the solutions of 
the former equation are not necessarily globally smooth for a positive e. A magnetic 
vortex, which reflects the behavior of the whole family (mg), is not the only type of 
singularity. Each function could independently nucleate the so-called harmonic 
bubbles. The question, whether the formation of bubbles really occurs in the Landau- 
Lifshitz equations, is open. In [10], the authors have to preclude the possible bubbling. 
Consequently, the conditions on the field v and on the degrees of the magnetic vortices 
are quite strict. We do not need to adopt these assumptions for O- We allow the 
degrees of both signs and we do not impose any restrictions on the structure of the 
external fields. Since we are not constrained by the physical considerations, the fields 
F and G are unrelated to each other. 

While treating the singularities of o, our overall strategy is the same as in the 
works m and m- We define the vortex positions as the concentration points of the 
Jacobian of u^- To identify the vortex motion law, we use the evolution equation for the 
Jacobian in combination with the estimates on the energy excess. Both major parts of 
the argument have to be adjusted to the situation when the equation contains forcing 
terms. 

The concentration properties of the Jacobian are tied to the Ginzburg-Landau energy 
of the solution. The behavior of the energy is the most complicated aspect of (|T|). The 
functional Egiug) can increase with time because the convective terms bring additional 
amount of energy into the system. This situation is new compared to the evolution 
Ginzburg-Landau equations. For them, the energy either decreases (heat and mixed 
flows) or is conserved (Schrodinger flow). We show that the energy can not increase 
too strongly. This requires a careful analysis of the impact that the convective terms 
have on the energy evolution. We find that the roles of the terms driven by F and 
G are quite different. To quantify the effect of the first of them, it is sufficient to 
apply the standard integral inequalities. For the second, the usual tricks do not provide 
satisfactory estimates. We rather make use of the precise structure of the emerging 
terms. 

We establish the vortex motion law by, roughly speaking, passing to the limit in the 
evolution equation for the Jacobian. The equation includes terms that depend on the 
external fields. We have to comprehend the limiting behavior of these new quantities. 
To this end, we show that the stress tensor and the momentum of concentrate at 
the vortex cores. These facts are not well-known, but they follow readily from the 
concentration properties of the Jacobian. Very similar results have been used in uni 
for the Landau-Lifshitz-Slonczewski equation. When we pass to the limit, we use, at a 
certain point, the estimates in terms of the energy excess. For this reason, we have to 
start with the well-prepared initial data. The energy of our initial data is approximately 
minimal among all functions with the given vortex configuration. 

In [^ and [53], equations similar to 0 appear in the context of the physically 
complete gauge-invariant Ginzburg-Landau theory. The equations stem from a special 
gauge transformation in a Ginzburg-Landau equation with magnetic field. The initial 
equation is, in addition, coupled with a boundary current. Exactly the boundary current 
produces the extra convective terms. The authors derive the effective motion law for 
these equations, but in regimes different from ours. In [23] . the vortex dynamics is 
mixed, like in Still, due to the scaling, it is driven by the external vector fields 
alone (no da^. W in the equation). The work [2tij treats the perturbed Ginzburg-Landau 
heat flow. In both cases, physics prescribes Neumann boundary conditions and some 
structural features of the external vector fields. In particular, the fields do not depend 
on time. 

The article is organized as follows. Section 2 contains all preliminary information on 
the notation and on the Ginzburg-Landau vortices. In the same section, we formulate our 
theorem and comment on the existence and regularity of the solution to ([T]). Section [3] 
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is dedicated to the analysis of the Ginzburg-Landau energy and of the kinetic energy of 
the solutions. In Section 01 we obtain the existence and continuity of the vortex paths. 
Moreover, we establish the results on the concentration. It remains to show that the 
vortex trajectories satisfy ©. We prove this via the Gronwall argument in Sections 
and [6] 


2. Preliminaries 

2.1. General notation. We distinguish between vectors in M" for n ^ 2 and scalars, 

either in R or in C ~ The vectors are in bold. Scalars and vector components are 
in light. Still, for u € C we write u = if we view u as a vector in 

For u,v € C we use the real scalar product {u,v) := + u^v^. The product of 

M G C with i is the rotation by tt/ 2: iu := [—v^,u^). The “vector product” for u,v G C 
is given hy u x v '■= {iu, v) = . 

The matrix Id is the identity matrix of size two. For two matrices A = (ajk) and 
B = {bjk), A: B is their Frobenius product A: B = Ujkbjk- 

The set Bs{x) is the closed ball of radius s centered at a; G R^. We denote by a an 
iV-tuple (oi, ...,ajsf) of distinct points belonging to the domain D. We set 

Pa := imn{i|afe - ai|,dist(afc,9II)}. 

The operator V is reserved for the differentiation with respect to the spatial variables. 
For a function i/i : —>■ R, curl i/i is a vector field curl i/! := {—d2'4’,diijj). For a vector field 

(p '■ D ^ R^, curie/) := is a scalar function. For a matrix A = (ajk(x)), 

div^l is the vector with the components divA^ = djUjk- For a scalar function ip, 
grad curl i/) is a 2 x 2 matrix 

grad curl V. = (■*;♦ ^ 

The space ((7°’“)* is the dual of the space of a-H61der continuous functions on D. 
The space is the dual of the space of a-Holder continuous functions that are 

compactly supported in D. The spaces (C'°’“([0,T] x D))* and (C'°’“([0,T] x D))* are 
defined in the similar way. Sometimes we shall need the elements of these spaces, and 
sometimes - vectors composed of such elements. We do not distinguish the two situations 
in our notation. 

We write = o(l) for a family of positive quantities q^, if lime_>o (?£ = 0. In the 
same manner, we write q^ ^ r^, if qe/ve = o(l). 

Throughout the text, C is a positive constant that does not depend on e. It may 
change from line to line. 

2.2. Quantities related to Ug. For any sufficiently smooth function 
u : R^ —>■ C, the Jacobian of u is given by 

(5) Jiu) ■= det(VM). 

The Jacobian is an exact 0-form. We have that J{u) = ^ curlwhere the vector 

(6) j{u) := {u X Vu) = (it X diu,u x d 2 u) 

is the current of it. In the same vein, for u = u{x, t) : R^ x R_|_ —>■ C, its total current is 
identified with the 1-form j^l = (it x djit) dt-|- (it x 9iit) dxi -I- (it x d 2 u) dx 2 - Accordingly, 
the total Jacobian of it is the exterior derivative of jti{u)- In other words, Jti{u) = 
dt A dxi +V’^dtA dx 2 + J{u) dxi A dx 2 , where J(it) is the usual Jacobian ([5]) and 
= dt{u X Uj) — dj{u X itt). We write {J{u), V{u)) for the three-dimensional vector 
representing Jti. The two-dimensional vector V (u) = is the velocity of the 

Jacobian. Since the form Jti is exact, J{u) and V{u) satisfy the equation 

(7) J(it) — curl V'(ii) = 0. 

The momentum of it is the vector p{u) := (9tit, Vit) = {{dtu,diu), {dtu,d 2 u)). 

The stress tensor is a matrix (Vit 0 Vit) with the entries (Vit 0 'S/u)jk = {djU, dku). 
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The energy density of u is the function 

Geiu) := + ^(1 - \u\^f. 

It is the integrand in the Ginzburg-Landau energy The rescaled energy density of 
the solution to (HD is given by 

( 8 ) g,e{t) ■■= keCeiUeit)). 


2.3. The vortices and the renormalized energy. The content of this section is 
standard and is taken from [I] and [S]. It mostly concerns the static case. 

For a function u : H —>■ C, a vortex is an isolated zero around which u/\u\ has a 
nontrivial winding number. The winding number is called degree. For e going to 0, the 
limiting vortices of the family (ug) are distinct points Ok & D with degrees dk G Z. The 
vortices of high degrees are believed to be unstable. We restrict our discussion to the 
case dk G {± 1 }. 

We complement our equation with either Dirichlet or zero Neumann boundary con¬ 
ditions. In the Dirichlet case 

(9) Ug\g^=g{x), 

the function g G C°°{dD,Ei^) has a nontrivial degree. This imposes a topological con¬ 
straint on the solutions of ©• A function u : D —> C that coincides with g on dD must 
have vortices in the interior oi D. In the Neumann case 

( 10 ) 

we assume that the initial data does have vortices. Moreover, we always assume that 
the initial data 


(11) Ue{x,0) = u°{x) 

is smooth. 

A vortex carries an unbounded amount of energy that is proportional to log(l/e). 
Therefore, the Ginzburg-Landau energy of a function with vortices diverges as s goes to 
0. The inequality 


(12) F;e(ue) ^ C'log(l/e) 

provides a non-sharp upper bound on the actual number of vortices. If (fT^ holds, the 
Jacobians of Ug are precompact in the strong topology of ((7°’“)* for every a > 0. They 
concentrate at the vortex sites [S]: 

N N 

J{Ug) TT^ dfcJaj, with TT^ Idfel < C. 
k=l k=l 


A vortex configuration (a, d) is the collection of vortex positions a together with the 
collection of corresponding degrees d G (±1)-^. We say that a family (ug) exhibits the 
vortex configuration {a,d), if J{ug) —^ '^'^k=i'^kSak- 

For a vortex configuration (a, d), the associated canonical harmonic map is an 

§^-valued map that solves the system 


(13) 


divj{u^.) = 0, 

curlj(it*) = 2TTj2k=idkSak- 


The boundary conditions are u* = 5 in the Dirichlet case and (t',j(M*)) = 0 in the 
Neumann case. In the former situation, the canonical harmonic map is defined uniquely, 
in the latter - up to a constant. 

The renormalized energy W of the vortex configuration (a, d) is given by 

O r ^ 

dx — d| log(l/s)y 

with It* = We shall sometimes write W{a) instead of W{a, d), since in our situa¬ 

tion the parameters d are fixed. In both Dirichlet and Neumann case, the renormalized 
energy has the form W{a) = ^kdj log \ak — aj \ + boundary terms. This shows 
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that the vortices interact like charged particles, with the degree playing the role of the 
charge. The renormalized energy blows up if two vortices bump into each other or a 
vortex hits the boundary of D. 

We denote by dakW{a) the vector ( 9^1 W(a), 9^2 W(a)). The identity below links the 
stress tensor of the canonical harmonic map to the gradient of the renormalized energy. 

Proposition 2.1 1 [TT| Proposition 3.2; [^, Lemma 8 ). Suppose that a smooth function 
if is affine in Bs(ak) for some 0 < s < p{a) and a fixed k € {1,..., N}. Suppose moreover 
that ip = 0 in a neighborhood of Oj for j ^ k. Let rt* be the canonical harmonic map 
associated to {a,d). Then, there holds 

(15) y*(gradcurlv?) : j{u^))dx = -cml(p{ak) • dai,W{a,d). 

D 

The Ginzburg-Landau energy and the renormalized energy are related as follows. 

Proposition 2.2. Suppose that {ug) exhibits a vortex configuration (a, d) with 
d G (±1)'^. Then there exists a constant 7 > 0 such that 

(16) liminf(£lE(uE) — 7rA^log(l/e)) > iVy + W{a, d). 

£->■0 

The generic constant 7 is the finite part of the energy of the vortex core. 

The approximate energy We{a,d) is given by 

Ws{a,d) := 7 rA^log(l/£) + iVy + W{a,d). 

It is the minimal Ginzburg-Landau energy that a family (ug) with a vortex configuration 
(a, d) can have. For such a family, the energy excess is defined by 

(17) D,:=E,{ue)-W,{a,d). 

If (me) depends on time, so does D^. 

The family (uf) is well-prepared, if for it = o(l). 

2.4. Requirements on the external vector fields. The vector fields F{x,t) and 
G{x,t) are smooth in both variables and satisfy 

= G{x,0)\g^ = 0 . 

Then, the compatibility conditions for the initial-boundary value problems O - ® - 
(ED and (E - CHD - (ED are fulfilled. In the Neumann case, we need in addition the 
equality 

(18) {G{x,t),i>)\QD = 0 

to hold for alH ^ 0. Here, iz is the unit outer normal on dD. 

2.5. Statement of the main result. The aim of the present work is to prove the 
theorem below. 


Theorem 2.1. Let (u^) be the family of solutions to CD with either Dirichlet or Neu¬ 
mann boundary conditions. Suppose that the initial data u® are well-prepared. In other 
words, there exists a configuration (a^,d) of N vortices with d G (±1)^ such that 

(19) = 7 riVlog(l/£) -I- Ny -|- W{a°, d) + o(l) 


and 

( 20 ) 


N 

'S dfeJjjO in the space ((7°’^)*. 

k=l 


Then, there exists a time T* > 0 such that for t G [0,T*), the following holds. 

(1) The family {ugit)) is well-prepared; 

(2) the Jacobian J{ue{t)) converges in (C)?’^)* to the measure 

(3) the vortex trajectories ak(t) solve the system of ODEs 


(Ao -\-i)ak = -^dakW{a{t),d) -\- F{ak{t),t) -j- dkiG{ak{t),t), 
afe(O) = al. 
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In the Dirichlet case, two vortices collide at time T*. In the Neumann case, at T* either 
two vortices collide or a vortex hits the boundary of the domain. 

2.6. On the existence of the solution to the equation ([T]). First, we can always 
find initial data satisfying the assumptions of Theorem l2.1l - see [5], Lemma 14. Second, 
the following proposition holds true. 

Proposition 2.3. For every fixed £ > 0, each of the problems (HJ - (H])- (fTTl) and © - 
(fTUl)- (fTT|) has a unique global classical solution. 

Proof. The reasoning is the same as in m, Remark 2.4. The equation is equivalent 
to a strictly parabolic system of semilinear equations on the real functions ul and 
For such a system, we construct a weak solution via the Galerkin approximation. We 
then improve the regularity with the help of the Calderon-Zygmund estimates m- □ 

Even though the system on ul and ul is parabolic, the weak maximum principle does 
not hold for it. This is due to the strong coupling. Hence, the uniform estimate |Me| ^ 1 
is not available, even in the Dirichlet case. The lack of this additional information is 
the only reason to use the excess energy estimates. Otherwise we could employ a more 
robust method of m for the derivation of the effective equation. This would allow us 
to lift the assumption on the well-preparedness of the initial data. 

2.7. Conservation laws. Here, we derive several intermediate evolution equations on 
the quantities defined in Section [531 

For the field {G ■ V)iUe, we easily find that 


(22) 

((G • V)iUe, Vrtg) = iG ■ J{ufi) 

and 


(23) 

{{G-V)tUe,dtu,)=-{G,Viue)) 


We obtain the conservation law for the energy density by multiplying both sides of 
o with dtUe and using (l23ll . The law is given by 

(24) dtCeiUe) = - ke{F,p{Ue)) + {G,V{Ue)) +divp{Ue). 

We also need the conservation law for the Jacobian. We first use the identity © to 
see that dtJ{ue) = cuTl{dtUe x Vu^) = cml (idtUe.'Vue) . We then substitute idtUe from 
©• Several formal transformations and (1551) finally yield 

(25) dtJ{ue) + AECurlp(u£) = curldiv (Vug 0 Vug) 

— fcgCurl {F ■ (Vue ® Vwg)) — curl {{iG) ■ J{ue)). 
The conservation law for the mass reads 

(26) + Ag(Mg X dtUe) + ke{F,j{ue)) + (G',V( ^~I""I )) = divj(Mg). 

We arrive at it when we multiply © with iug. 

3. Energy estimates and corollaries 

In this section, we study the Ginzburg-Landau energy and the kinetic energy of the 
solution to ©• We first prove a rather crude estimate on the Ginzburg-Landau energy 
in Lemma l3.II From it, we deduce an optimal estimate on the kinetic energy in Propo¬ 
sition o This implies the relative compactness of the total Jacobian in Gorollarv l3.ll 
We finally improve the control on the Ginzburg-Landau energy in Proposition 13.21 

Lemma 3.1. For every T > 0, there exists a constant C > 0 such that for allO ^ t ^ T, 
there holds 

(27) EMt)) ^Clogil/e). 

Proof. We integrate (153) over the set D x [0,T]. Then, we obtain 

T T T 

(28) E,{u,{T)) + X, JI \dtu,\^ = E,{ul) - k, I J{F,p{u,)) + J J {G,V{u,)). 

0 D 0 D 0 D 
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We estimate the integrals on the right-hand side of (051) . To the first of them, we apply 
successively the Cauchy-Schwartz inequality and the Young inequality with appropri¬ 
ately chosen coefficients. This gives 


(29) 


ke 



1 

2 


T T 

-\-Cke J Ee{Ue{t))dt. 

0 D 0 


The estimate for the second integral is slightly more involved. The decisive step is to 
use the structure of V{ue) first. The remaining estimates are again standard. We start 
with integration by parts 


T T T 

J J {G,V{Ue)) = -^ J J divG ■ {Ue X dtUe) + ^ J J {dtG,j{Ue)) 

0 D 0 D 0 D 

+ ^J (G(0), j(u°)) - {G{T),j{u,{T)))dx. 

D 

The integral over dD, which appears when we integrate by parts with respect to x, is 
equal to zero. In the Dirichlet case, this follows from the fact that the boundary data 
g{x) is independent of time. In the Neumann case, this follows from the additional 
assumption (ITSll on the field G. 

We work with the three new terms. We take into account the bound 

T T T 

(30) J Jk\^ = J j (Kix, t)\'^-1)+TC\D) + Ce^ j E,{u,{t))dt. 

0 D 0 D 0 

Here, is the Lebesgue measure on We obtain, by using the Cauchy-Schwartz 
inequality, o, and the Young inequality, the estimates 


div G ■ {Ue X dtUe) 


0 D 


T T 

J +Ce^\og{l/e)- J E^{ueit))dt + C\og{l/e), 

0 D 0 


T 


0 D 


{dtG,j{Ue)) 



1 

2 




J E^{Ue{t)) dt + Cs‘^ log(l/e) J Ee{ue{t)) dt + Clog(l/^), 
0 0 


and 


j{G{Q),j{ul)) - {G{T),j{u,{T)))dx 

D 

^ C\og{lle) + log(l/£) + ke)E,{ue{T)). 
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Together with (1^ . this gives 


(31) EMT)) + f I I \dtu,\^ ^ E,{ul) + C{e^ log(l/e) + k,)E,{u,{T)) 


0 D 


1 

+ C{e^ log(l/e) + ke) J Es{ue{t)) dt + C'log(l/£). 


By choosing e small enough, we absorb the term C{e^ log(l/e) + ks)Es{ue{T)) on the 
left-hand side. Now, the Gronwall Lemma and the well-preparedness of the initial data 
imply (l?7l) . □ 

Proposition 3.1. There exists a time T > 0 and a constant C > 0 such that 

T 

(32) 


ke J j \dtUe\^ dxdt ^ C. 


0 D 

Proof. Our proof is by contradiction. It is a modification of the argument in Lemma 
III.l of [21]. We assume that the claim is false. Then we rescale time in such a way that 
the vortices do not move, but the energy is dissipated strongly enough. This contradicts 
the well-preparedness of the initial data m- 

Suppose that the statement of the proposition does not hold. Then there exist se¬ 
quences —>■ 0 and 0 such that 

in 

ke^ / / 

0 D 

For the rest of this proof, we write instead of Ue„ and kn instead of . 

We first rescale time by kn and define the functions Vn{x, t) := Un{x, knt). For them, 
our assumption gives the equality 


(33) 


\dtVn? = 1, 


0 D 


with T„ ;= t„log(l/e„). We rescale time once again by t„ and define the functions 
Wn{x,t) ■= Vn{x,Tnt). Note that Wn{x,t) = Un{x,tnt). Foi' Wn, there holds 


(34) 

and 


i 

J J \dtWn\‘^ =Tn^ l0g(l/e„) 


0 D 


(35) Ee{Wn{t)) = Ee„{Un{tnt)) ^ Clog(l/£„). 

Due to Theorem 3 in [22], the estimates dMj) and ([351) imply that the total Ja¬ 
cobian of Wn is relatively compact. Hence, the limits lim„_>.oo =/(rCn(t)) = Jft) and 
lim„^oo V'(r(;„(t)) = V(t) exist in ((7°’“)*. For them, we have dtJ{t) — curiy(t) = 0 
for every t G [0,1]. In other words, the vorticity of Wn is transported by V. Moreover, 
V satisfies the product estimate from m 


(36) 


1 


jj.v.x 

0 D 


^ lim inf 

n—foo 


1 


0 D 


\X • VWn 



1 

2 


The inequality holds for every Holder continuous scalar function ip and every Holder 
continuous vector field X. The right-hand side of (|551) is zero, by (IMl) and (1551) . Con¬ 
sequently, V = 0, and the vortices of Wn do not move on the time interval [0,1]. They 
are the same as the vortices of and J{wn{t)) for every t G [0,1]. 
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We return to Vn- Its vortices do not move until the time t„. With (ITBl) . we have 
Es„{vn{Tn)) ^ 7rA^log(l/£) + Nj + W , d) + o(l) = E!;{u°) + o(l). Moreover, 

(37) Es„{Ve{Te)) = Es{Ue{te)) ^ Clog(l/e), 

by Lemma l3.II 

The functions solve an equation of the same form as O- The only difference is 
the scaling of the time-derivative. Therefore, as in the proof of the Lemma [23 we find 
that 


(38) Ao y j \dtVn\^ = Ee^{u°^) - Ee^{Vn{Tn)) 

0 D 

Tn Tn 

-kr^J J{F{k^t),p{Vr,))+J J{G{knt),V{Vn)). 

0 D 0 D 


The quantity on the right-hand side in the first line is either positive of order o(l) or 
negative. Below we show that both integrals in the second line are of order o(l). This 
leads to a contradiction, which concludes the proof. 

For the first of the integrals we use the Cauchy-Schwartz inequality, ((33ll and (l37ll to 
obtain 


k 


n 


{F{knt),p{Vn)) 


< Ckn 



0 D 


0 


o(l). 


In the second integral, we switch to the time-scale of and get 

Ffi 1 

J J {G{knt),V{Vn)) = J J{G{tr,t),V{w„)). 

0 D 0 D 


As we already know, V{wn{t)) —?> 0 for t G [0,1]. Then, this integral is of order o(l). 

□ 


The most important consequence of the energy estimates established so far is 


Corollary 3.1. The total Jacobian {J{ue),V{ue)) of is relatively compact in the 
space ((^“’“([O,T] x D))*, for every a G (0,1]. 

Proof. The claim follows directly from (EH) and (12^ by the virtue of Theorem 5.1 in 
[5]. Even though the Theorem is formulated for the space ((^^’“([O, T] x D))*, the proof 
works for the smaller space (C°’“([0,T] x D))* as well. See a related discussion in the 
proof of Theorem 3 in [22]. □ 


Proposition 3.2. There exists a constant C > 0 such that, for all t G [0, T], there holds 
(39) Es{usit)) 7rA^log(l/£) -I- C. 


Proof. We return to ([22)1 and set T = T. Now we can estimate the integrals containing 
F and G much more precisely. With ([27| and (Ea), we obtain that 


ke 





For the second integral, we use Corollary 13.II Since V{ue) is relatively compact and 
G is smooth, there exists a constant C such that 

T 


(40) 


J J{G,V{u,)) 


CC. 


0 D 


This gives the inequality Es{ue{T)) ^ E^{vP^) + C. Now, (|3^ is actually a consequence 
of ([11]). □ 
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4. Vortex trajectories and concentration 


The vortices are the concentration points of the Jacobian. From Corollary 13.11 we 
deduce in Proposition id. li the existence and continuity of the vortex paths. We study the 
concentration properties of the stress tensor and the momentum of Ug in Propositions 
HU and 14.51 Proposition 14.61 is a technical fact, which we need in the next section. It 
follows from the energy bounds. 

Proposition 4.1. Let T be the time found in Provosition \S.1\. For every sequence 
e„ —>■ 0, the sequence {J(ue^),V{ug^)) converges, modulo a subsequence, in the space 
T] X D))*. For the limit {J{t),V{t)), we have 

N N 

J{i) = andV{t) = -ttj ^ 

fc=l k=l 

The functions fk(t) ■ [OiT] —>■ belong to the space F[^{0,T; D). The coefficients dk 

are integers, with dk G {±1}. 

Proof. The existence of a convergent subsequence is given by Corollarv lJ.il The measure 
J{t) is, for every t, a weighted sum of deltas by Theorem 3.1 of [5]. Moreover, Theorem 3 
of [22] yields that J{t) belongs to the space T], ((7°^)*). The regularity of the 

functions fk is determined in m, Propositions III.l and III.2. The last two facts also 
imply that the coefficients dk are constant. Finally, the form of V(t) is completely 
specified by the representation for J{t) via ([7|). □ 


In addition, we have a mobility bound on the functions fkit). 


Proposition 4.2 (\2‘2\. Corollary 7 ). Let fkit) be the functions found in Proposition 
\4.1\ Then, for all 0 ^ ti < t 2 ^ f, there holds 


(41) 



dt < liminf kg 
6 —^0 


//ia.T. 


ti D 


We use the part of Propositionthat concerns V (ug) in a different, but equivalent, 
form. Namely, we have that for all 0 ^ < ^2 ^ f and all smooth vector fields 

w : D X [0, T] —>■ R.^, there holds 


(42) 


lim 
6 —^0 


*2 ^ t2 

j J {y{u,{t)),w) = TT^dk j {ikit),iw{ik{t),f))dt. 

ti D ti 


The set of paths ^{f) coincides with the solution set a{f) of (1211) at the time t = 0. 
This is due to the condition (1^0)1 on the initial data. Our aim is to prove that the two 
sets coincide for all t £ [0,T*). We point out that ^{t) is produced by some sequence 
En —>■ 0. Up to the end of Proposition 16.11 we use the notation e —>■ 0 for this particular 
sequence. Eventually, Proposition 16.11 and Theorem o imply that the family ^{t) is 
independent of the sequence £„ —>■ 0. 

In view of gu, we may assume that T is smaller than the collision time T*. Then, 
by Proposition 13.21 and the definition of the energy excess (ITT)) , we find that 


(43) De{t) < log(l/e) 

for every t G [0, T]. 

We now deduce from Proposition 14.II the remaining results on the concentration. 


Proposition 4.3 (Concentration of the rescaled energy density). Suppose that J{ue{t)) 
converges as described in Proposition \4 b\ Then the rescaled energy density /ig (<) ([8|) 
concentrates at the vortex cores: 

N 

(44) lim^e(t) = 7r^%(t) m ((7°’^)*, 

/c=l 


for every t G [0, T]. 
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N 

Proof. See [3], Theorem 1.4.4. We denote the limiting measure tt 

fc=i 


Proposition 4.4 (Concentration of the stress tensor). For every t G [0, T] and for every 
smooth vector field v : I? —>■ there holds 


(45) 


N 

lim / V ■ kg{'S/Ue{t) G^Vue{t)) dx = v{^kit)). 

6^0 j 

D 


Proof. Due to ProDOsition l3.21 the family {ke{yue{t)(SiVug{t))) is relatively compact in 
the space (C^d)*. Moreover, due to Proposition |331 for every Borel set A C Il\ {^(t)}, 
there holds 

J ke\{djUe{t),diUe{t))\ dx < J dx = o(l). 

A A 

Therefore, the stress tensor converges to the measure of the form 
where Ak(t) are some 2x2 matrices. It remains to show that Ak(t) = 7rld, for all k 
and t. To this end, we use the result of Kurzke and Spirn on the equipartition of the 
Ginzburg-Landau energy (Proposition Ed])- 

We fix t, set a = and apply Proposition lA.ll in for a fc S {1, ...,N}. 

The assumptions are satisfied for all sufficiently small e, due to Propositions 14.II and l4. 31 
Moreover, we may take as Kq the constant C from Proposition 13.21 Consequently, the 
constant is the same for all t and k. We finally obtain 


/ Ak{t) — irld dx ^ lim / ke{Vue (Si\7ue) — irlddx 

£^■0 J 


^ hm — , 


= 0 . 


□ 


Proposition 4.5 (Concentration of the momentum). For all 0 ^ ti < t 2 ^ T and all 
smooth vector fields w : D x [0,T] ^ there holds 


*2 


(46) 


lim 

S —^0 


{keP{Ue{t)),w) = -TT 


ti D 


N 

E 


{ikit),w{^k{t),t))dt. 


Proof. The family {k^p{us)) is precompact in the space {C°’^{D x [0,T]))*, grace to 
(I32]) and (|39p . Moreover, due to (l44ll . the limiting measure-valued vector field p has the 
form p = TT Y.k=i Pk{t)dik{t)- The functions p^. : [0, f] R^ satisfy /J |pfc| • |^fe| dt < oo 
for every /c G {1,..., N}. It remains to prove that Pf. = —TT^k almost everywhere. 

We observe that for any smooth function (j){x,t), there holds 


(47) 


t2 

- J J Vcj) dp = J (I){x,t2)dfi{t2) 

ti D D 


J (j){x,ti)dfi{ti) 

D 


t2 


dt(j){x^ t) dfi{t) dt. 


ti D 


Indeed, we multiply (|2^ with A:e0(x,t) and integrate over D x [^ 1 ,^ 2 ]- We obtain 


t2 t2 t2 

\,k, j j \dtu,\^Tj> + kl j J{F,p{uSm-keJ J {G,V{u,))cjT 

ti D ti D t\ D 


= -k 


t2 t2 

J J dtee{ue)cl)-\- ke J J divp('Ue)0. 


D 


D 


All terms in the first line converge to zero as e goes to 0. This follows from (15^ . (I^Hl) . 
and (gOl)- We integrate by parts in the second line and arrive at (HZl). 

We construct a special test function (j){x,t)- For po ■= inf{p^(t),t G [0,T]} > 0, we 
set cr = po/2. We fix a smooth function (p(x) such that <p(x) = x in 73^(0) and ip{x) = 0 
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outside of B2cr(0). We finally set 4‘{x,t) := for a fixed k S A^}. 

When we consider (H71) with the test functions (fP, j = 1,2, the identity takes the form 

-/ PUt)dt = Tr f 

Jti Jti 

The equality holds for all [^ 1 ,^ 2 ] C [0,T]. Thus p^. = —ir^k almost everywhere, as 
required. □ 

Proposition 4.6. The quantity drvj{ue) converges to 0 in the dual space 0/x 

[ 0 ,T]). 

Proof. We test the conservation law for the mass (ES with a function (p from Wq^’^(D x 
[ 0,T]) and integrate over D x This gives 


T T T 

-J J = J J (p ■ + ke J j (p-[F,j(Ue)) 

0 D 0 D 0 D 

T T 

+ 1J + J J ct>-{u,xdtu,). 

0 D 0 D 

We estimate the terms on the right-hand side using Propositions 13.11 and 13.21 together 
with the Sobolev embedding Wq’'^{D x [0,T’]) L°°{D x [OjT]). We have 


11 

0 D 


11 

0 D 


^ C||5t0||i2ev'log(l/e); 


1 1 
k,J J cP{F,jiu,))\^CkM\L^{ J J 


0 D 


0 D 



0 D 





i 1 

J J p{G,vC-^))\ = \J |v(</>G)-(i^))| <C||V<^||^c4ev'bi(lM 


0 D 


0 D 


1 

J J P{Ue X dtUe) ^ XeCWPWh^ ijj ^ ' (/ / 


0 D 


0 D 


0 D 




All the quantities on the right-hand side are of the order o(l), and so the claim follows. 

□ 


5. Estimates involving the energy excess 

We measure the difference between a{t) and ^{t) on the level of the energy in terms 
of the energy excess. We now introduce two more quantities that we use for the same 
purpose. 

The first of them is the vortex position error ri{t). It is a vector in with the 
two-dimensional components rikft) = Cfe(t) — ak(t). The second quantity is the ODE 
error R{t) G We define it via the components as 

(48) Rk{t) := Xoik{t) + dkiik{t) P \da^W{a{t)) - F{ak{t)) - dkiG{ak{t)). 

In other words, Rkit) = (Aq -I- dki)fik{t). Consequently, 

(49) \R{t)\ = ^Xl + l\rj{t)l 
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In this section, we establish a relation between Dg{t), ri{t) and R{t). It is given by the 
inequalities (ISTl) and (15^ . 

We first need the following auxiliary estimate. 

Proposition 5.1. For every 0 ^ ti < t 2 ^ T, the energy excess satisfies the inequality 

N 

(50) DS 2 ) - DSi) ^ 

k=l 

Flere, the functions Q^it) := F(£^klt),t) + dkiG{fk{t),t) are vector-valued and the func¬ 
tions Skit) := /q VP(afe(t)+s-77fe(t), <)+(ifcV(iG)(afc(t)+s-77/c(t), t) ds are matrix-valued. 

Proof. By the definition of the energy excess, we have 

Deit 2 ) - Dsiti) = Es{us{t 2 )) - Esiueih)) + W{a{ti)) - W{a{t 2 )). 

We rewrite the right-hand side with the help of (PT|) and (PSI) . We obtain that 

De{t2) — Dg{ti) = Bi -\- B2 -i- B3. 


*2 

J -2Ao(dfc, ?)fe) + iQk,Vk) + {Sk ■ dk, dk) dt -\- o(l). 


The terms Bj are given by 


B^=- 


*2 JY *2 

Ae J J\dtUe\'^-\- ttXoY J\dk\‘^dt, 

ti D 


*2 Jkf *2 

B 2 = -ke J j iF,p{ue)) - ttY j iFiak{t),t),dk)dt, 

ti D 


and 


B. = 


*2 pf *2 

dk J{iG{akit),t),dk)dt. 

t. n fc=l i. 


tl D 


We consider them one by one. 
Proposition 14.21 implies that 


Bi ^ ttAo ^ / {\ak\^ - dt + 0 ( 1 ) = ttAo^ / {-rjk^Vk + 2dfc) dt + o(l) 

tl tl 

N ‘2 

<-27rAo^ idk,rik)dt -\- o{l). 

k —11 

tl 

The estimates for B 2 and B 3 follow from similar arguments. By Proposition 14.51 we 
obtain 


B 2 =t^Y [i^i^k{t),t),ik) - iF{akit),t),dk)dt-\-oil) 

k —11 

tl 

N 

= ^Y J W’ dk) + (F(Cfe(i), <) - Fiakit), t),dk) dt + o(l). 


*1 


We represent the vector F(^fe) — F{ak) with the help of the Mean Value Theorem by 
F{^k{t),t) - F{ak{t),t) = (/q^ VF{ak{t) + s ■ r]k{t),t)ds) ■ r]k{t). 
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For i? 3 , we use and the Mean Value Theorem to get 


Bs = n'^dk ({iG{^k{t),t),f]k)dt 

N 

+ 7r^dfc / {jQV{iG){ak{t) + s-rikit),t)ds)-rik(t),dk)dt + o{l). 

fc=i i 

We sum the estimates for Bj and arrive at (IHTIl) . □ 

Proposition 5.2. For every T < T, there exists a constant C = C{T,aP) such that, 
for 0 ^ ti < t 2 ^ T, there holds 

*2 *2 

(51) Dsit 2 ) + \'n{t 2 )\ < Ds{ti) + |r 7 (ti)| +C J \r]\dt + C j |i2| dt + o(l). 

ti ti 

Proof. We estimate the right-hand side of (1501) in terms of R. The time T is strictly 
smaller than the collision time for the vortex motion law m- Therefore, the right- 
hand side of m is bounded on [0,T]. The functions Qf. and Sk are bounded on this 
interval as well, because the fields F and G are smooth. Hence, there exists a constant 
C = {T,aP) such that 

sup |d(t)| ^ C, max sup \Qk{t)\ ^ C, and max sup |S'fc(t)| ^ C. 

tG[0,T] ^ te[0,T] ie[0,T] 

With this C and 0 ^ ti < t 2 ^ T, the inequality dsni) simplifies to 

*2 *2 

De{t 2 ) De{ti) + + C j \'n(t)\dt + o(l). 

ti ti 

We add to both sides of the inequality the term \r]{t 2 )\ — \r]{ti)\. On the right-hand 
side, we estimate this term from above. We have ||T 7 (t 2 )| — |’ 7 (ii)|| ^ 

On the other hand, we have that |^|’ 7 (t)|| = |’ 7 (t)|. This way, we see that 

t2 t2 

De{t 2 ) + \'n[t 2 )\ - |’7(ii)l ^ Defti) + C j \v{t)\dt + C J \r]{t)\dt + o(l). 


This inequality together with yields (|5T|) . 


□ 


The proposition below is the last major ingredient of the proof of Theorem 12.II The 
estimate has been proven in m with the tools developed in The argument relies 
on the energy estimates, the concentration of the Jacobians, and the fact that divj(ue) 
converges to zero. Grace to Propositions 13.2114.11 and 14.61 the proof works in our case 
as well. Therefore, we do not repeat it. 

Proposition 5.3 (|11]. Proposition 5.4 ). There exists a constant K > 1 depending only 
on D and N with the following property. Suppose that at a time to ^ 0, D^ito) = o(l) 
and T](to) = 0. Let t > 0 be so small that kg fjj \dtUg\‘^ ^ C. Set a := ^Pa(to)- 
Then for every S ^ t such that Dg(t) ^ C and |i7(t)| ^ cr/(4iF) on the interval [to, to+<5], 
and for every (p € Wo'°°{D) that is affine in each Ba{ak{to)), there holds 


(52) 


tQ-\-6 


to D 


(gradcurl(/?) : ((Vug (g) Vue) — (j(w*) j(“*))) dx 

to-\-S 


dt 


^ C y Ds{t) + \r]{t)\dt + o{l) 


, = u^:{a{t),d). 


with U: 
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6. Proof of Theorem 12.II 


In this section, we finish the proof of Theorem 12.11 We first establish a short-time 
result in Proposition 16.11 With its help, we finally obtain the full statement on the 
time-interval [0,T*). 

Proposition 6.1. There exists a time r > 0 such that the conclusions of Theorem \2.1\ 
hold on [0, r]. 


Proof. The time T is strictly smaller than the time of the first collision for the solution 
to (ED- Therefore, the number A^, := sup{|d(t)|,< G [0,T]} is well-defined. We set 
p* := s-nd a := p*/2. The radius a is strictly positive. For 

the constant K > 1 from Proposition 15.111 we define tq := min{T, cr/(4111 A*)}. Then, 
ak{t) G i?CT/( 4 R')(a°), for every k G {1,..., N} and all t G [0, tq]. Since the paths are 
continuous, there exists a time ri > 0 such that ^k{t) G for every k G {1, ...,iV} 

and alH G [0, ri]. We finally set 


T min{To,Ti}. 

We show that Theorem 12.11 holds true on the time interval [0 ,t]. 

We fix a cut-off function x<yi^) ^ such that xi^) = 1 in .6^(0) and xi^) = 0 

outside of B2a{0)- We fix a number k G {l,...,iV} and define a vector-valued test 
function cj) = {cj)^,cj)^) with 4>{x) = (ix) ■ Xcr{^ ~ In the ball B^{a'j.), there holds 

4>^{x) = —X 2 , (t>^{x) = xi, curl(}()^(x) = (1,0), curl(/)^(x) = (0,1). 

Outside of the ball i32cr(a)!), all these functions are equal to zero. 

We apply Propositions 14.II and 14.51 to the test functions cjA, j = 1,2. We obtain 


dt / J(ue(t))(jA dxdt - \e / {cml 42 ,p{ue)) dxdt ^ TT / dt{dk{i£.ky 


Xoil) dt. 


D 


0 D 


For the equation error 

+ dkiii + -da,,W{a{t)y - F^{ak{t),t) - dkifGy {ak{t),f), 

TT 

there holds 


r T 

(53) TT J lirn J dlt J J{ue{t))4y dx — J {curl 4>y p{ue)) dx 


0 D D 

+ dakW{a{t)y - TTF^[ak{t),t) - ndkiiCy {ak{t),t) dt. 


We substitute the first two terms on the right-hand side of (IM)) from (1^ . This yields 
the inequality 

r 

TT J |ii^| dt < lnn(di + I2 + I3), 

0 

where the terms Ii, I 2 , I 3 are given by 


h = 


42 ■ curldiv {Vue®Vue)dx + da^W {a{t)y 


dt, 


0 D 



0 D 


I3 = 


{cnrlfy ,iG ■ J(ue)) dx — ndkiiGy {ak{t),t) 


dt. 


0 D 

We estimate them one by one. 
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By the choice of cj), da^W{a{t)y = curlc/)-^ • da^W{a{t)) in B„{a^). Therefore, by 
Propositions 12.11 and 15.31 we have 


h 


(gradcurl) : ((Vite 0 Vm^) - (j(M*(a)) 0 j(u*(a)))) 


0 D 


^cj Dyt) + \v{t)\ + o{i). 
0 


For the term I 2 we obtain, by Proposition 14.41 and the smoothness of F, that 


I2=7T J - F^{akit),t) + 0 ( 1 ) J \r]k\ + o(l). 

0 0 

Analogously, by Proposition 14.II we have that 

r T 

/a = TT y* {iGy{^k{t),t)-{iGy{ak{t),t) dt F o{l) ^ C J \rjk\dt F o{l). 
0 0 
We sum over j and k to see that 

r r 

J |i?| dt ^ linr C J \ri{t)\ + 0^(1) dt. 


We combine this inequality with (l5T|) . At t = 0, 77 ( 0 ) = 0 and Dg(0) = o(l). With the 
Gronwall Lemma we conclude that r]{t) = 0 and L?e(t) = o(l) for all t G [0, r]. The first 
identity establishes the vortex motion law. The second shows that the solution remains 
well-prepared on [0,t]. □ 


Proof of Theorom, \2. 1[ completed. We iterate our argument starting at the time r, which 
is found in Proposition 16.II We apply successively Propositions 13.II - 16.II on a new time- 
interval [r, T-ff] with f > 0. It remains to check that we can reach, by iterating so, the 
terminal time T*. 

We argue by contradiction. Let T be the supremum over all times for which The¬ 
orem [51T] holds. By Proposition 16.11 T is positive. Suppose that T < T*. We shall 
show that we can actually apply Proposition 16.II at the time T. This will contradict the 
assumption on maximality of T. 

We have to prove that fkiT) = ak{T) and Ds{T) = o(l). For every T G [0,T), (1^ 
gives 

T T T 

(54) j j\dtu,y = E^ul) - E^uyT)) - k, j j{F,p{uy)+J J (G^Viu,)). 

0 D 0 D 0 D 

Using Propositions 16.li H31 and 14.11 we see that the right-hand side is equal to 


W{a^)-W{a{T)) 



— TT 


N 

/ {iG{ak(t),t) 

fc=l n 


, d/c) -l-o(l). 


The paths Ok remain distinct for T = T because T < T*. Thus, the quantity above 
remains bounded. This means that /J< C . By Theorem 3 of [H], the 
functions fk are continuous on the closed interval [0,T]. Since a(T) is equal to ^(T) on 
[0,T), and both collections of functions are continuous, we obtain a(T) = ^(T). 

We now consider (1541) on the interval [T, T] instead of [0, T], Estimating the integrals 
with F and G as before, we get 


E^uF)) < FFT)) 


— TT 



{F{ak{t), f) + dkiG{ak{t), t) 


,dk) + o(l). 
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By the well-preparedness of (ue) at all times T <T, we have that 

E,{ue{T)) = log(l/e) + TVy + W{a{T)) + o(l). 

For each k G we have by (I4ip and the regularity of F and G that 

f ft 

I j{F{ak) + dkiG{au),au)\^C{ j \G{ak)\^ F \F{ak)\^)^ ■ { j\ak\^)^ = C {f - T)\. 

T T T 

This gives 

E,{ue{f)) \og(l/e) + N-1 + W{a{T)) + C{f-T)\+o{l). 

Since W{a{T)) —>• W{a{T)) for T —^ T, we obtain that 

lim D^{f) ^ lim lim {W{a(T)) + Cif - T)\ F o(l) - W{a{f))) = 0. 

£—>•0 £—>-0 T^T 

The proof is now complete. □ 


Appendix A. Equipartition of the Ginzburg-Landau energy 


Proposition A.l ([E], Theorem 1 and Proposition 15). Suppose that for some 
a > 0, we have 





and 

Then there holds 
on Kq. 


j e^{Ue) 


dx < Trlog —\- Kq. 
e 


f ke(Vus^Vue) — 7rlddx ^ with KI depending continuously 
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